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Yan events to probe Z-boson production mechanisms. 
In quantum chromodynamics (QCD) at order a s this 
occurs either through the annihilation process with a 
gluon (G) in the final state (qq — > j*/Z G), or via 
the Compton process with a quark in the final state 
(qG — > 7*/Z q). The emission of final state q/G gives 
7*/Z transverse momentum [l[ (we define the produc- 
tion P T = Pt{i* /Z) = P T (e + e~) before final state radi- 
ation). 

The general expression for the angular distribution Q 
is described by the polar (8) and azimuthal (<f)) angles 
of the decay-electron in the Collins-Soper (CS) frame |3|- 
When integrated over cos# or <fi, respectively, the decay- 
electron angular distribution is described by: 

da oc (1 + cos 2 0) + -A (l ~ 3 cos 2 ff) + A 4 cos 6 (1) 
a cos 9 2 



da 



cx 1 + fa cos cf) + (3% cos 2(f) + 07 sin <j) + /3 5 sin 2<j) (2) 



where p 3 = 3ttA 3 /16, (3 2 = A 2 /4, (3 7 = 3ttA 7 /16, /3 5 = 
Ac, /A. The Aq and A4 are extracted from Eq. [IJ and A2 
and A3 are extracted from Eq. [2j while A5 and A7 are 
expected to be zero 0. 

Perturbative QCD (pQCD) makes definite predictions 
for the angular coefficients A), 2, 3, 4 (A) and A2 are the 
same for 7* or Z exchange, and A3 and A4 originate from 
the j* /Z interference). For the qq — > 7* jZ G annihila- 
tion process pQCD at order a s predicts that the angular 
coefficients Aq and A2 are equal 0-0] and can be ana- 
lytically described by Af = Af = P*/(M 2 e+e _ + P 2 ) 
(Eq. 3). At higher order, there are small deviations from 
the above expression (Eq. 3) which depend on PDFs and 
dilepton rapidity (y) 

For the qG — > 7* jZ q Compton process, Aq and A2 
depend on parton distribution functions (PDFs) and y. 
However, in pQCD at order a s , when averaged over y, 
Aq and A2 are approximately described [a la] by Aq = 
Af a 5P 2 /(M 2 +e _ + 5P 2 ) (Eq. 4). 

At order a s , the Lam- Tung relation (A = A 2 ) [l(| is 
valid for both qq and qG processes [jj. Fixed-order pQCD 
calculations at order a 2 [2| as well as QCD resummation 
calculations to all orders 6], indicate that violations of 
the Lam- Tung relation are small. The Lam- Tung relation 
is only valid for vector (spin-1) gluons. It is badly broken 
for scalar (spin-0) gluons . Therefore, confirmation of 
the Lam- Tung relation is a fundamental test of the vector 
gluon nature of QCD and is equivalent to a measurement 
of the spin of the gluon. A previous determination of 
the gluon spin was made from a study of 3-jet events 



qq G) in e + e annihilation [12j. 



To date, the Lam- Tung relation has been tested only at 
fixed-target experiments using samples of low mass Drell- 
Yan dilepton pairs at relatively low transverse momen- 
tum. In this region, non-perturbative higher-twist effects 
can be significant [l3L Il4j . Some experiments report large 



violations @, 14- III, and one experiment [l?} is consis- 
tent with the Lam- Tung relation. Here we report on the 
first test of the Lam- Tung relation at large dilepton mass 
and high transverse momentum, where non-perturbative 
higher-twist effects are expected to be negligible. 

Fixed order pQCD calculations Q and Monte Carlo 
(MC) simulations at next-to-leading order (NLO) (e.g. 
dyra d fl8| and madgraph [19], and pythia in Z+ljct 
mode [20(| ) indicate that there is a significant (w 30%) 
contribution of the Compton process to the production of 
7*/Z bosons at the Tevatron. Therefore, as shown in Fig. 
[21 these calculations yield values of A and A 2 which are 
larger than the pure annihilation process prediction (Eq. 
3). Similar results are predicted by POWHEG [2l|, a NLO 
MC with additional parton showering, and fewz [22I ] 
which is a next-to-next-to-leading order (NNLO) QCD 
calculation. 

In contrast, the default, LO version of pythia [23|, 
and VBP [24[ (an MC generator based on QCD resum- 
mation) predict values of Aq and A2 which are close to 
Eq. 3 (which is only correct if the qq process is dom- 
inant). The RESBOS [25( MC generator, which is also 
based on QCD resummation, predicts values of Aq and 
A 2 close to Eq. 3 at low Pt, and larger values (close 
to the predictions of fixed order pQCD) at high Pr, as 
shown in Fig. [3J Therefore, measurements of Aq and A2 
as a function of Pt elucidate the relative contributions 
between the annihilation and Compton processes. 

In this Letter, we report on the first measurement of 
the angular coefficients Aq, A2, A3 and A4, for pp — > 
7*/Z — > e + e~ + X events in the Z boson mass region 
(66 < M ee < 116 GeV/c 2 ) produced at ^/s = 1.96 TeV. 
We also report on the first test of the Lam- Tung relation 
at high transverse momentum. 

The sample used corresponds to an integrated lumi- 
nosity of 2.1 fb _1 collected by the CDF II Detector at 
Fermilab ^ during 2004-2007. Charged particle direc- 
tions and momenta are measured by an open-cell drift 
chamber (COT), a silicon vertex detector (SVX), and 
an intermediate silicon layer in a 1.4 T magnetic field. 
Projective-tower-geometry calorimeters and outer muon 
detectors enclose the magnetic tracking volume. The cov- 
erage of COT tracking in pseudorapidity is \rj\ < 1.2 [l[. 
Reconstructed tracks are used to determine the pp col- 
lision point along the beam line, which is required to 
be within z = ±60 cm of the center of the detector. 
The energies and directions [l[ of electrons, photons, and 
jets are measured by two separate calorimeters: central 
(\rj\ < 1.1) and plug (1.1 < \r]\ < 3.6). Each calorime- 
ter has an electromagnetic compartment with a shower 
maximum detector followed by a hadronic compartment. 
Three topologies of e + e~ pairs are considered: two cen- 
tral electrons (CC), one central and one plug electron 
(CP), and two plug electrons (PP). Events with at least 
one electron with high Et are selected online. Offline 
refined selection requires the electron to have Et > 25 
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FIG. 1: Di-electron Pt spectrum of data, default (CDF 
Tuned) pythia prediction, and backgrounds (QCD and elec- 
troweak process). The mass range corresponds to 66 < M ee < 
116 GeV/c 2 . 



GeV for CC and PP events, and E T > 20 GeV for CP 
events in the fiducial regions of the calorimeters, the cen- 
tral (|j7„| < 1.1) and plug (1.2 < \r) e \ < 2.8). To minimize 
background, the second electron candidate is required to 
have E T > 15 GeV for CC, E T > 25 GeV for PP, and 
Et > 20 GeV for CP events. The selection criteria listed 
above are the same as in the related previous publication 
27j of the Z rapidity distribution, but are augmented in 
this analysis with the additional requirement that both 
electrons have an associated track in the SVX. The data 
sample consists of about 140 000 events. The fractional 
contribution of the total QCD background (2-jet events 
misidentified as a Drell-Yan pairs) to the number of se- 
lected events is 0.3%. This is determined by studying the 
distribution of transverse energy in a cone surrounding 
the center of the electromagnetic cluster in the calorime- 
ter. The total background from electroweak (WW, WZ, 



MK+jets, and Z — > r + r ) and tt processes is estimated 
from simulation to be 0.2%. 

The effect of the acceptance on the angular distribu- 



tions is modeled using the PYTHIA MC generator [23 1 
combined with a GEANT [28[ simulation of the CDF de- 
tector. The pythia generator includes a LO QCD in- 
teraction (qq j*/Z), initial state QCD radiation, par- 
ton shower fragmentation, the 7* jZ — > e + e~ decay, and 
photon radiation from the final state. The version of 



pythia used at CDF has additional ad-hoc tuning 23] 
(referred to as default pythia) in order to accurately rep- 
resent the 7* jZ boson transverse momentum distribution 
measured in data. Further tuning was introduced in or- 
der to ensure that the MC simulation correctly described 
the rapidity, as well as the correlations between rapid- 
ity and transverse momentum that are observed in the 
data. To reconstruct the simulated events in the same 
way as data, the calorimeter energy scale, resolutions, 
and selection efficiencies used in the detector simulation 
are tuned [27| using data. Figure Q] shows the di-electron 




cose,. 



FIG. 2: The cos6» distribution of data and default (CDF 
Tuned) pythia prediction. 



Pt spectrum for data, the default pythia prediction, 
and the backgrounds. There is good agreement between 
data and PYTHIA prediction. Figure [5] shows the cos 9 
distribution for data and the default pythia prediction 
and its ratio. 

The analysis is performed in five bins of transverse 
momentum as shown in Table |U For each transverse mo- 
mentum range, data and MC simulated events are binned 
in cos 8 and <f>. The MC events are re- weighted to gener- 
ate the expected angular distributions (cos 9 and <p) for 
a range of values of Aq and A4, and Ai and A3, respec- 
tively. The angular distributions from the re-weighted 
MC events are compared to the data in the reconstructed 
level and the angular coefficients which give a maximum 
log-likelihood value are determined as the best coeffi- 
cients to describe the data. The Aq an d A4 are deter- 
mined by the comparison of the data to MC distributions 
in cos 6 and the Ai and A3 are determined in The 
normalization factor of the data to MC events is also in- 
cluded as one of fit parameters. The results are shown in 
Fig. [3j and in Table Q] with statistical and systematic un- 
certainties. The correlation between extracted values of 
Aq and A2, A3 and A4 is negligible. The systematic un- 
certainties originating from backgrounds, electron identi- 
fication efficiency, SVX tracking efficiency, boson Pt and 
rapidity modeling, and modeling of detector material are 
considered. The dominant source is the background esti- 
mate. Most of systematic uncertainties are discussed in 
reference 27[ and the effect of these uncertainties on the 
shape of the angular distribution is small. 

The data are in good agreement with the Lam- Tung 
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FIG. 3: Comparison of the measured values of Ao, A2, A3 and 
An (for 66 < M ee < 116 GeV/c 2 ), shown with statistical and 
systematic uncertainties combined in quadrature, to theory 
predictions. The data are plotted at the mean Pt of the events 
for each bin. The last bin corresponds to Pt > 55 GeV/c 
with no upper limit. The horizontal uncertainty is RMS of the 
transverse momenta in each bin. Agreement [29( is found with 
the predictions of fewz and powheg (shown) , and also with 
dyrad , MADGRAPH , and pythia Z +l-jet MC (not shown). 
The data do not favor [2!j the predictions of default pythia , 
and vbp . Also shown are the pure qq — > 7* /Z G annihilation 
diagram prediction and the qG — > J*/Z q Compton process 
prediction from the pythia Z +l-jet MC. 



relation Aq — A2 = 0, which is expected in QCD with 
vector gluons. The values of A — A 2 for the five Pt 
bins are 0.00 ± 0.03, 0.04 ± 0.05, 0.03 ± 0.07, 0.02 ±0.11, 
and 0.01 ± 0.14 (statistical and systematic uncertainties 
combined), which average to (Ao — A 2 )=0.02 ± 0.02. At 
low Pt the measured values of A and A 2 are well de- 
scribed by the qq — ► 7* jZ G annihilation function (Eq. 
3). At high Pt the larger values show that both the 
annihilation and Compton processes contribute to the 
cross section [2(|. Our results are in agreement [29| 
with fixed-order perturbation theory calculations includ- 
ing DYRA D [lij . MADGRAPH [lj|, PYTHIA Z+l jet [2p] |. 

POWHEG |2l|. and fewz [22| (all of these give similar 
predictions). We find that the values of A 3 and A4 are 
in agreement with the predictions of all models (A4 is 
calculated with sin 2 6*w = 0.232). 

In summary, we present the first measurement of the 
angular coefficients in the production of 7*/Z bosons 
at large transverse momenta, and the first test of the 
Lam- Tung relation at high transverse momentum. We 
find good agreement with the predictions of QCD fixed- 
order perturbation theory, and with the Lam- Tung re- 
lation Aq = A 2 . The measurements presented here are 
statistically limited. An analysis with larger samples in 
both muon and electron channels is currently under way. 
A comparison of these results with future measurements 
at the LHC would provide additional tests of production 
mechanisms since the contribution of the Compton pro- 
cess (qG — >• 7* jZ q) at the LHC is expected to be larger. 



TABLE I: The measured angular coefficients (measured value 
± stat. error ± syst. error). The mean Pt of the events 
in the five bins are 4.8, 14.1, 26.0, 42.9, and 73.7 GeV/c, 
respectively. 



Pt bin 


A) (xHT 1 ) 


A 2 (xlO" 1 ) 


0-10 
10-20 
20-35 
35-55 
> 55 


0.17 ±0.14 ±0.07 
0.42 ±0.25 ±0.07 
0.86 ±0.39 ±0.08 
3.11 ±0.59 ±0.10 
4.97 ±0.61 ±0.10 


0.16 ±0.26 ±0.06 
-0.01 ±0.35 ±0.16 
0.52 ±0.51 ±0.29 
2.88 ±0.84 ±0.19 
4.83 ± 1.24 ±0.02 


Pt bin 


M (xlO- 1 ) 


A 4 (xl0 _i ) 


0-10 
10-20 
20-35 
35-55 
> 55 


-0.04 ±0.12 ±0.01 
0.18 ±0.16 ±0.01 
0.14 ±0.24 ±0.01 
-0.19 ±0.41 ±0.04 
-0.47 ±0.56 ±0.02 


1.10 ±0.10 ±0.01 
1.01 ±0.17 ±0.01 
1.56 ±0.26 ±0.01 
0.52 ± 0.42 ± 0.03 
0.85 ± 0.50 ± 0.05 



We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda- 



7 



tion; the A. P. Sloan Foundation; the Bundesministerium 
fur Bildung und Forschung, Germany; the Korean World 
Class University Program, the National Research Foun- 
dation of Korea; the Science and Technology Facilities 
Council and the Royal Society, UK; the Institut National 
de Physique Nucleaire et Physique des Particules/CNRS; 
the Russian Foundation for Basic Research; the Ministe- 
rio de Ciencia e Innovation, and Programa Consolider- 
Ingenio 2010, Spain; the Slovak R&D Agency; and the 
Academy of Finland. 



[1] In the CDF detector frame, the positive z axis is de- 
fine by the proton beam direction. The pseudorapidity 
is rj — — lntan(#/2). For an e + e~ pair Pt = Psind, 
Et = EsinO, where 9 is the polar angle between the 
particle direction and the z axis. For an e + e~ pair, 
y — | In g^p* , where P and P z are the magnitude and 
z-component of the momentum, and E and M e + e - are 
the energy and mass. 

[2] E. Mirkes and J. Ohnemus, Phys. Rev. D 50, 5692 (1994); 
Phys. Rev. D 51, 4891 (1995). 

[3] J. C. Collins and D. E. Soper, Phys. Rev. D 16, 2219 
(1977). The CS frame is the center of mass of the "y* /Z, 
where the z axis is defined as the bisector of the p and 
p beams. The polar angle is defined as the angle of the 
positive lepton with respect to the proton direction. 

[4] J. C. Collins, Phys. Rev. Lett. 42, 291 (1979). 

[5] D. Boer and W. Vogelsang, Phys. Rev. D 74, 014004 
(2006). 

[6] E. L. Berger, J. Qiu, and R. A. Rodriguez-Pedraza, Phys. 

Lett. B 656, 74 (2007). 
[7] A. Bodek, Eur. Phys. J. C 67, 321 (2010). 
[8] S. Falciano et al. (NA10 Collaboration) Z. Phys. C 31, 

513 (1986). 

[9] J. Lindfors, Physica Scripta 20, 19 (1979). 
[10] C. S. Lam and W. K. Tung, Phys. Lett. B 80, 228 (1979). 
[11] Ao — A2 ~ —2 in the Gottfried- Jackson frame, see N. 

Arteaga-Romero, N. Niclaidis, and J. Silva, Phys. Rev. 

Lett. 52, 172 (1984). 
[12] R. Brandelik et al. (TASSO Collaboration) Phys. Lett. 

B 97, 453 (1980) (their data was consistent with spin 1 

gluons, but 3.8 standard deviations from the prediction 

for spin gluons). 
[13] D. Boer, Phys. Rev. D 60, 014012 (1999). 
[14] J. S. Conway et al. (E615 Collaboration), Phys. Rev. D 

39, 92 (1989). 

[15] M. Guanziroli et al. (NA10 Collaboration) Z. Phys. C 37, 



545 (1988). 

[16] J. G. Heinrich et a/.(E615 Collaboration), Phys. Rev. D 
44, 1909 (1991). 

[17] L. Y. Zhu et al. (FNAL-E866/NuSea Collaboration), 
Phys. Rev. Lett. 99, 082301 (2007). 

[18] W.T. Giele, E.W.N. Glover, and D. A. Kosower, Nucl. 
Phys. B 403, 633 (1993). 

[19] F. Maltoni and T. Stelzer, J. High Energy Phys. 302, 27 
(2003). We use CTEQ6.6 PDFs for madgraph . 

[20] T. Sjostrand et al, J. High Energy Phys. 05, 26 (2006) 
We use pythia 6.4 Z+l jet matrix element event genera- 
tor (MSEL = 13) with parton showering (MSTP(61) = 2 
and MSTJ(41) = 1) and CTEQ5L PDFs. We also gen- 
erate pythia 6.4 Z + 1 jet without any parton showeing 
(MSTP(61) = and MSTJ(41) = 0). The pythia Z + 1 
jet without parton showering gives the same angular co- 
efficients with pythia Z + l jet with parton showering. 

[21] We have used the most recent version of POWHEG. 
See S. Alioli, P. Nason, C. Oleari, and E. Re, 
larXiv:1009.5594V l (2010), with CTEQ6.6 PDFs. 

[22] K. Melnikov and F. Petriello, Phys. Rev. Lett. 96 231803 
(2006); Phys. Rev. D 74 114017 (2006), with MSTW2008 
NNLO PDFs. 

[23] Default (CDF tuned) pythia (6.216) includes the qq and 
part of qG (t-channel only) processes.. T. Sj0strand et 
al., J. High Energy Phys. 05, 26 (2006). We use the 
default (MSEL=11) LO matrix element (Z + jet) 
with CTEQ5L PDFs. The parton showering produces 
the boson P T . The CDF EWK/TOP standard W/Z P T 
tuning parameters are: MSTP(91) = 1, PARP(91)=2.10, 
PARP(93) = 15 for the low Pt Gaussian smearing, with 
PART(62) = 1.25 and PARP(62)=0.2 for the P T evolution 
in 7-25 GeV region. The underlying event is included as 
Tune A. The QED parton showering uses the same ma- 
chinery as QCD parton showering aside from coupling 
differences. 

[24] R. K. Ellis and S. Veseli, Nucl. Phys. B 511, 649 (1998). 

We use CTEQ6.6 PDFs with vbp . 
[25] F. Landry, R. Brock, P.M. Nadolsky, and C.-P. Yuan, 

Phys. Rev. D 67 073016 (2003). We use CTEQ6.6 PDFs 

with RESBOS . 

[26] D. Acosta et al. (CDF Collaboration), Phys. Rev. D 71, 
032001 (2005). 

[27] T. Aaltonen et al. (CDF Collaboration) Phys. Lett. B 
692, 232 (2010); Jiyeon Han, Ph.D. Thesis, University 
of Rochester, FERMILAB-THESIS-2008-65. 

[28] GEANT: CERN Program Library Long Writeup W5013. 

[29] For the four highest Pt bins the \ 2 values for Ao are 11.7, 
9.4, 5.6 and 3.7 for default pythia (qq and t channel qG), 
vbp, RESBOS, and fewz (qq+qG), respectively. 



